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ABSTRACT 

In order to use strong gravitational lens time delays to measure precise and accurate cosmological 
parameters the effects of mass along the line of sight must be taken into account. We present a 
method to achieve this by constraining the probability distribution function of the effective line of 
sight convergence K ox t- The method is based on matching the observed overdensity in the weighted 
number of galaxies to that found in mock catalogs with K cxt obtained by ray-tracing through structure 
formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, 
and redshift. This additional information reduces the uncertainty of K oxt from a K ^0.06 to ~0.04 for 
very overdense lines of sight like that of the system B1608+656. For more common lines of sight, 
a K is reduced to <0.03, corresponding to an uncertainty of < 3% on distance. This uncertainty has 
comparable effects on cosmological parameters to that arising from the mass model of the deflector and 
its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to 
constrain K oxt almost as well as with spectroscopic redshifts. As an illustration, wc apply our method 
to the system B1608+656. Our most reliable K ox t estimator gives ct k =0.047 down from 0.065 using 
only galaxy counts. Although deeper multi-band observations of the field of B1608+656 are necessary 
to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to 
characterize the immediate environment, is an effective way to increase the precision of time-delay 
cosmography. 

Subject headings: gravitational lensing — methods: numerical — (cosmology:) cosmological parame- 
ters 



1. INTRODUCTION 

Strong gravitational lensing can be used to study 
a number of important astrophysical topics, rang- 
ing from the mass distribution within galaxies (e.g., 
IKochanekl 119911 : iTreu fc Koox )mansl I2002D and clusters 
of galaxies (e.g.. [Kneib et al.lll993t ILimousin et al.|[2u07l : 
1 12008T) to magnifying distant sources (e.) 

hard et alBOQS) 



Sand et al.1 izuusii to magmiymg Q istc 
Pettini et al.ll2000t [Ellis~et al.ll2001b iRic 



to cosmography (e.g., IKochanekl 119961 : iSchechter et al 



. apny (e 

1997 ; [Qgurill2007D (see recent reviews by iCourbin et al.l 
20021 iMevlan et al.[|2006[ : ITreu! 120101: iBartelmannl 120101: 
Kneib fc Nataraianll201ll for extensive references) 



The power of lensing as a tool for astrophysics stems 
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from its ability to measure mass independent of its dy- 
namical state, as well as from its ability to magnify 
background sources. Detailed mass models of the de- 
flector galaxies, groups, or clusters of galaxies can be 
determined by reproducing the strong lensing observ- 
ables, i.e. multiple ima ge time delays, positions, and 
fluxes (e.g.. iKc cton 201(| and references therein). How- 
ever, one of the main limitations of gravitational lens- 
ing is due to th e so-called "mass-sheet degeneracy" 
(IFalco et all fl985t ISchneider fcSeTO fl995t fSaTial [27)0(1 
IWucknitzH2002Tr This emerges from considering the solu- 
tion to the lens equation, characterized by a certain sur- 
face mass density k for the deflector in units of the crit- 
ical density S C rit- By linear transformations, one finds 
that an infinite family of solutions can be obtained. The 
family of solutions results in a range of inferred prop- 
erties both for the mass model of the main deflector as 
well as for the properties of the lensed source. One needs 
to break the mass-sheet degeneracy using physical argu- 
ments in order to fully exploit the power of gravitational 
lensing. 

A number of strategies have been adopted to break 
this degenera cy both in the strong and weak lensing 
regimes (e.g.. iBroadhurst et al.lll995 ; iKneib et aLl 120031 : 
iBradac et all l2004t iSonnenfeld et ail 1201 ID A common 
strategy is requiring the surface mass density of the main 
deflector to vanish at large distances from its center. 
This solution is physically equivalent to assuming that 
the distribution of mass along the line of sight (LOS), 
excluding the plane of the main deflector, is uniform and 
equal to the average of the Universe. This approxima- 
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tion is sufficient for many applications of gravitational 
lcnsing, resulting typically in uncertainties of only a few 
percent in the inferred mass distribution of the main de- 
flector and in the luminosity and size of the lensed source 
(iSeliakJ 119941: IKeeton et al.lll997t iKoopmans et all 120061: 
iTreu et all 120091 : iHoekstra et al.ll2011[ ). For higher pre- 
cision measurements however, one needs to determine 
the effects of the line of sight mass distribution. It is 
customary to condense these effects into an equivalent 
additional mass sheet at the redshift of the main de- 
flector with uniform surface mass density K e xt, which 
can be positive or negative depending on whether the 
line of sight is over or u nderdense with respect to the 
average of the Universe ()Schneiderl Il997f ). In practice, 
mass-sheet degeneracy can be broken, i.e., the value of 
K ox t can be inferred, by constraining in an independent 
way (1) the mass of the lens galaxy via stellar kine- 
matics fe.g.JTreu fc Koopma ns 2002: lKoopmans fc Treu 



matics fe.g.JTreu fc Koopmansll2002t iKoopmar 
20031 : iBarnabe et al.ll2009HAuger et alj |2010t IS 
2010: ISonnenfeld et all 120121) . and/or (2) 



iuvu et al 



2010t ISonnenfeld et all 120121 ). and/or (2) the total 
mass of any intervening mass structures along the 
line of sight via imaging an d spectroscopy of objects 
along the line of sight fe.g.. IKeeton fc Zabludofil 120041: 
Fassnacht et al.ll2006t iMomcheva et alJl200e " ' 



20101: iWong et al.ll20nnFassnacht et alJl201lD ~ 
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Determining re ext is especially important for doing 
precision cosmography with gravitational lens time de- 
lays. Gravitational lens time delays are the differ- 
ence in the arrival time of photons along the paths 
corresponding to multiple images arising from ge- 
ometric and general relativistic effects. For vari- 
able sources, like active galactic nuclei, time delays 
can b e measured via dedica te d monitoring campaigns 
(e.g.. IFassnacht et al I 120021 iFohlmeister et all 120081: 
Paraficz et al.1 120091: iCourbin et al.l 120111: iTewes et al.l 



20 VI'.. If the mass distribution in the plane of the de 



flector and along the line of sight is known, the mea- 
sured time delays can be converted to the so-called time- 
delay distance D&t, a combination of three angular di- 
ameter distances (e.g., ITreu! [2010h . In turn, the time- 
delay distance contains information on cosmolo gical pa- 
rame ters, primarily the Hubble constant Hq ( Refsdall 
1 19641). but also the curvature and other parameter s 
( Coe fc Mousta kas 2001 ISuvu et alJl2010t ILinderll20lll) . 



To first approximation, time delays constrain cosmology 
in a similar way to baryon acoustic oscillation experi- 
ments and therefore are an excellent complement to other 
probes like the cosmic microwave background a nd type 
la Supernovae (|Linderll20Hl: iDas fc Linderll2012l ). It has 
been shown that with current techniques and sufficient 
ancillary data a single gravitational lens m easures the 
time-delay distance to approximately 5-6% (jSuvu et al.1 

l2TnTil2T)i3l ). 

In cases like B1608+656 where the time delays are 
known to ~2% and the mass model of the main deflec- 
tor and its immediate environment is exquisitely con- 
strained by the data, the dominant source of uncertainty 
is the distribut ion of mass along the line of sight (e.g., 
iBar-Kanall 19961 ) , i.e. effectively the external convergence 
K ext . Specifically, -Da* oc (1 — K C xt) _1 - Thus, for small 
ftextj the uncertainty a K translates directly into relative 
uncertainty in time-delay distance, currently dominating 
the overall error budget. 



ISuvu et "ail (|2010i ) constrained K ox t by using the den- 
sity of galaxies within 45" of B1608+656, which was mea- 
sured to be twi ce that of an ave r age fie ld observed at the 
same depth bv IFassnacht et ahl (|2011| ). Then, by select- 
ing lines of sight with the same galaxy over density Co-ai in 
the M illennium ray-tracing simulations of iHilbert et al.1 
(|2009|) . they measured the conditional probability dis- 
tribution function (PDF), -P(/c e xt|Cgai), which was then 
used as a prior in deriving cosmological information. Sev- 
eral efforts are underway to explore ways to further re- 
duce this source of uncertainty. Those include spectro- 
scopic and pho tomet ric surveys (pMomcheva et al.l 120061 : 
iWilliams et all 120081: lAuger et all 120071) as well as mea- 
surements of the external shear (|Suvu et al.ll2013| ) 

In this paper we focus on refining and developing prac- 
tical tools to estimate external convergence by comparing 
the output of cosmological numerical simulations with 
readily available observables such as those that can be 
derived from a galaxy photometric catalog. We define 
the relative overdensity £ as the value of selected observ- 
ables related to the stellar mass, the redshift, and the 
projected distance on the sky (angular distance) divided 
by the mean of the same o bservables over all lines of 
sight. Extending the work of ISuvu et all (pOlOl) we con- 
sider a set of observables with relative overdensities £ in 
addition to £gai, the number of galaxies along a LOS di- 
vided by the average, seeking observables that minimize 
<t k . In a companion paper, Collett et al. (2013, sub- 
mitted) a halo-model approach is used to perform a full 
reconstruction of the mass distribution along the line of 
sight, and thus provide a theoretical counterpart to the 
methods and weighting schemes developed here. 

The main result of this paper is that using these statis- 
tics it is possible to reduce significantly a K , using a 
multiband photometric catalog. The amount of gain de- 
pends on the specific line of sight, but for a system like 
B1608+656 it is possible to reduce it from 6 to 5% (even 
4% in the best cases). We note that even reducing the 
uncertainty on a single lens from 5 to 4% is extremely 
useful given how rare these strong lenses are and how 
time-consuming it is to obtain ancillary data like time 
delays and high-resolution images. Oversimplifying for 
the purpose of illustration, assuming a target precision 
of 1% from a sample of lenses, only 16 systems would 
be needed if each were precise to 4%, as opposed to 25 
systems at 5% precision eacl0. 

The paper is organized as follows. In Section [2] we 
briefly summarize the Millennium Simulation that forms 
the backbone of our study. In Section [3] we revisit the 
galaxy count statistic Cgai- In Section 2] wc introduce 
other statistics involving stellar mass, luminosity, red- 
shift, and distance. In Section [5] we test the influence of 
using photometric redshifts instead of spectroscopic red- 
shifts in our method. In Section [6] we apply our method 
for several statistics to the gravitational lens B1608+656. 
Throughout the paper magnitudes are given on the Vega 
scale. 

2. SUMMARY OF THE MILLENNIUM SIMULATION 

11 In reality, since lenses and sources will have different red- 
shifts, the likelihoods of the cosmological parameters from each lens 
will not be perfectly aligned in multiple dimensions, and therefore 
the combined uncertainty should decrease faster than 1 /yjV (e.g. 
Coc & Moustakas 2009; Dobke et al. 2009; Linder 2011). 
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Numerical simulations of large-scale structure provide 
a way to determine statistically the amount of external 
convergence for a lens syste m. By ray trac i ng th rough 
the Millennium Simulation (jSpringel et al.l I2005T ) , one 
can produce a simulated image of the sky populated with 
galaxies and a corresponding map of the external con- 
vergence, Kext, fo r a given source redshift ([Hilbert et al.l 
2007, 2008, 2009). The distribution of K ex t values asso- 
ciated with lines of sight in the simulation that resemble 
the lens system of interest (e.g., in terms of the over- 
density of galaxies around the lens system) provides a 
statistical estimate of the K ex t for the lens system. 

We use 64 simulated fields of 4 x 4 deg 2 from the Millen- 
nium Simulation containing galaxies with positions, red- 
shifts, magnitudes (in SDSS u,g,r,i,z and 2MASS J,H,K) 
an d stellar masses from the semianalytic galaxy model 
of IGuo et al.l ^OlllP^l Each field has approximately 5 
million galaxies with redshifts up to ~ 3.2 and an asso- 
ciated map of the external convergence on a 4096 x 4096 
grid for a source redshift of 1.4, typical of strong lensing 
systems like B1608+656. We use each position on the 
K ext map as a line of sight, and therefore have approxi- 
mately 10 9 lines of sight for the 64 fields. iHilbert et al.l 
(|2009f) and lSuvu etall (|2010f ) showed that the distribu- 
tion of Kext from strong lens lines of sight are very similar 
to the distribution for all lines of sightF 5 !. Therefore, we 
consider all lines of sight from the 64 fields in our analy- 
sis. These lines of sight provide the pool from which we 
select subsets that have observational properties (based 
on the galaxies) matching those of the lens system for 
determining the At oxt of the lens. 

3. GALAXY NUMBER COUNTS AS A PROBE OF K cxt 

In this section we revisit th e galax y overdensity statis- 
tic introduced bv lSuvu et al.l (|2010[) . We recall that the 
use of relative overdensities - instead of absolute densi- 
ties - in both the data and simulations is meant to min- 
imize the impact of theoretical and observational sys- 
tematic uncertainties. In particular this should reduce 
sensi tivity to the cho ice of a specific reference simula- 
tion (jSuvu et al.ll2010t ). 

In practice, we study -P(/v e xt|Cgai)i where £ ga i is given 
by the number of galaxies within 45", iV ga i, in a specific 
line of sight, divided by the average value for all lines of 
sight, iVgai; i.e. Cgai = For the Millennium Simu- 

lation data, N sa \ is readily computed by multiplying the 
total number of galaxies by nr 2 /A where r = 45" and 
A is the angular area of the entire field. However, in 
practice one must be careful of masked regions and edge 
effects. The choice of the radius is motivat e d by p ractical 
reasons. As discussed by IFassnacht et al.1 (|2011| ). this is 
typically the maximum radius that can be surveyed for 
a target in middle of one of the two chips, i.e. the stan- 
dard pointing of the Wide Field Camera of the Advanced 
Camera for Surveys (ACS) on board the Hubble Space 
Telescope (HST). We also impose restrictions similar to 

12 Obtained from the Millennium Database 

(Lemson & Virgo Consortium 2006). 

13 The Koxt constructed from the Millennium Simulation ex- 
cludes the convergence at the primary lens plane (i.e., the redshift 
of the lens galaxy) since this is already accounted for in the lens 
galaxy mass modeling. Therefore, the external convergence is truly 
external to the lens and is due to line-of-sight contributions. 



observations of B1608+656 such that all galaxies must 
have < z < z sourco (w here z SOUYCe = 1-394 is the source 
redshif t of B1608+656 , IFassnacht et all 119961). and fol- 
lowing EssnachFer^I]jl2^^^aTOT8^ < magi < 24.5. 
The results found by ISuvu et al.l (|2010D used only the 
latter constraint. In Section O we revisit this field with 
redshift information. Once the PDFs of K ox t are com- 
puted we define the width of the PDF a K as the semid- 
ifference of the 84 and 16 percentiles of P(re e xt|Cgai)- We 
remind the reader that in this paper we use the con- 
vergence maps detailed in Section [2] to supply our K cxt 
values. 

The K ex t PDF for a LOS with a known relative overden- 
sity in galaxy count is readily computed by first counting 
iVgai for each line of sight on the convergence maps. The 
convergence maps are defined on a regular grid with res- 
olution of ~3.5", and the grid points thus conveniently 
serve as locations of the lines of sight for our galaxy num- 
ber counts. Then, it is sufficient to select lines of sight 
with overdensity close to the desired value Cgai- In prac- 
tice, we select all lines of sight with 



Kal-CgalM <E (1) 

where E, the interval width, is some integer value. As we 
increase E, the number of lines of sight satisfying Equa- 
tion ([!]) also increases. As a general rule the parameter 
E should be chosen to be as small as possible while still 
leaving sufficient lines of sight to generate the PDF. Note 
that sample variance noise is already introduced by the 
simulations so we expect that varying E while keeping it 
smaller than ~ y/N gSL \ should not affect our results. 

A number of subtleties must be taken into account 
when constructing the « ex t PDFs. As expected, the dis- 
tributions are in general highly skewed. For example, 
fewer lines of sight become available at higher relative 
overdensities (and similarly at lower undcrdensities; for 
conciseness we shall discuss explicitly only the overden- 
sities in our examples). Because the number of lines 
of sight A^los a t a given galaxy count iV ga ] arc always 
greater than that at N ga \ + 1, simply constructing a PDF 
from all K ox t whose corresponding N Ra i satisfy Equation 
{J) would be biased towards lower JV ga i, and their respec- 
tive Kext values. We define K™ t d as the median K oxt for a 
given PDF, and find that a good indicator of the bias in 
a sample is the change in k™ t d as function of the interval 
width E. The narrowest interval E = 1 gives the closest 
estimate to the true median. Since lower N ga \ pull the en- 
tire PDF away from that with a galaxy count of Cgai-^gaij 
we expect that a notable change in K™ t d will occur when 
we increase E. This is shown in the top panel of Figure 
[TJ To circumvent this shift in K™® d that would lead to a 
bias in the time-delay distance determination, we weight 
each value of n cxt by N ^ for its respective iV ga i - that 
is, we find -/Vlos for the N ga x responsible for contributing 
a particular ft e xti and multiply by the inverse. There- 
fore each of the 2E n e xt PDFs for a given iV ga i within 

CgaiAgai ± E carries equal weight, and re ex t from higher 
iVg^ offset those from lower, thus ensuring that our over- 
all distribution remains relatively static. This is indeed 
the case, as the bottom panel of Figure [T] shows At™ j d 
declines much more slowly as a function of uncertainty 
than in the upper panel. 
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E (Galaxy Counts) 

Fig. 1. — Change in K™° d (solid) and a K (dashed) with E for 
Cgal = 0.75, 1.25, 1.5, and 2 as measured from numerical simula- 
tions. The change is computed for both (a) unweighted and (b) 
weighted PDF combination, by subtracting re™° d and <r K for the 
given distribution at E by those at E = 1. Unweighted refers to 
simply averaging K cx t values for all LOS that satisfy (JJ, whereas 
weighted indicates weighting inversely by the number of LOS of a 
particular galaxy count. The K™j d line for f ga i = 0.75 in panel (b) 
follows a nearly identical line to that of 1.25 and therefore is barely 
visible. 




Fig. 2.— PDFs for K ex t constructed from lines of sight that 
satisfy Equation JTJ for f ga i = 0.75, 1.25, 1.5, and 2 and E = 2 as 
measured from numerical simulations. 



However, maintaining a steady median has resulted in 
increasing the width of the distribution a K . This is ex- 
pected because PDFs on either side of CgaiA ga i are cen- 
tered slightly above or below our target. For this reason 
- and that of residual effects on the median - we see that 
lower galaxy interval widths (E) offer the best results, as- 
suming they can encompass sufficient data to adequately 
reduce statistical uncertainty. For the purposes of this 
paper unless otherwise noted we set E = 2. 

As an illustration of our method, in this paper we inves- 
tigate the undcrdcnse case £ ga i = 0.75 and the overdense 
cases Cg a i = f .25, f .5, and 2 (used in Figure [T]), although 
our method carries over to any arbitrary value of £ ga i. 
We show in Figure [2] the PDF of K ex t for these £gai val- 



ues. As expected, larger £ ga i produces a shift in the PDF 
along positive K C xt, however, it also increases a K . In other 
words more overdense lines of sight have higher conver- 
gence but also a broader range o f poss ible convergences. 
In the study of iFassnacht et all (|2011h . BI608+656 pos- 
sessed a relative overdensity of 2.18 without redshift cuts. 
We will thus look at the case of £gai = 2 with particular 
interest. 

4. ALTERNATE CHARACTERISTICS AS WEIGHTS 

Defining £ ga i = A ga i/A ga i is useful for constructing 
PDFs of K cxt for strong lenses. Counting iV ga i for any 
line of sight is straightforward and requires only that 
each object within 45" have z < Source and a flux greater 
than an observational limit. However, by using -/V ga i we 
neglect characteristics of an object that may play a sig- 
nificant role in gravitational lensing (i.e. mass, redshift, 
angular offset). By using quantities that are closely re- 
lated to lensing we expect that we should be able to 
reduce the uncertainty on K e xt- For example, we do not 
suspect all LOS with £ ga i = 2 to have exactly the same 
physical characteristics as B1608+656, hence we can use 
relative overdensities in observable features besides A ga i 
to construct an even tighter - and perhaps more relevant 
- PDF. 

4.1. Redefining £ for New Characteristics 

To identify LOS with overdensities with particular fea- 
tures, we need to design a weighting scheme such that all 
objects are not equivalent, but weighted by a feature rel- 
evant to lensing. By summing each object multiplied by 
its weight within 45" we define the weighted sum 

i=l 

where qi is the weight for object i. Note that for our def- 
inition of -/Vgai, qi — 1 for all i and therefore the weighted 
number of galaxies W ga .i simply corresponds to N ga \. We 
next define ( q = W q /W q as before. In accordance with 
Cgal- we wm set Q q = 0.75, 1.25, 1.5, and 2 for the purpose 
of demonstrating our general method. 

4.1.1. Weighting by Radius 

The angular separation between the source and a 
nearby object is a significant factor in distorting, and 
thereby shaping, the path through which the source's 
light passes. We expect then that each galaxy within 
45" does not contribute equally, but that those nearer 
the optical axis of the lens are more influential than 
those farther away. In particular for isothermal total 
mass distributions (appropriate for massive galaxies or 
around the scale radiu s of halos, e.g.. lGavazzi et al]l2007l : 
lLagattuta et al.ll2010|) , we expect the convergence to de- 
cline as the inverse projected distance from the deflec- 
tor. This gives rise to our first weighting method of 1/r. 
We will scale all objects that satisfy 10" < r < 45" by 
1/r. At r < 10" weighting becomes sensitive to small 
changes in r. Thus for r < 10" we allow each object to 
carry a weight of 1/10, giving us a continuous weight- 
ing function. We note that in general objects that are 
very close to the main deflector are more likely to be 
physically associated with it and exert a stronger impact 
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Fig. 3. — Contour plots showing f ga j vs. f 9 for various weighting methods, as measured from numerical simulations. Dark to lighter 
shades refer to 1,2,3, and 4<r regions, respectively. The dotted line marks £ ga i = ( q . We see that the different £ 9 are correlated, i.e. 
under/overdense LOSs in one metric are typically under/overdense in other metrics as well. However, they are not identical, implying that 
imposing multiple £ g conditions adds information and thus sharpens the PDF. 

though their halos do. For M n with n > 1 we will con- 



(e.g-. lKeeton fc Zabludofill2004D . Therefore it is prudent 
to model them explicitly, rather than considering them as 
part of the statistical lines of sight effect. This might re- 
quire obtaining as m uch information as poss i ble on them, 
inclu ding redshifts (jMomcheva et all 120061 : lAuger et al.l 
I2007T ) and possibly stellar velocity dispersions, especially 
if they are consistent with massive galaxies. 

4.1.2. Weighting by Redshift 

Objects close to the source along the line of sight have 
minimal lensing effects from the scaled deflection of light 
rays. Likewise, those nearest to the observer are rel- 
atively insignificant. In order to approximately account 
for this, our next heuristic weighting method is quadratic 



in z and defined as z s 



z 2 where 



and z are 



the redshifts of the source and the object along the LOS, 
respectively. For simplicity the notation of this weighting 



4.1.3. Weighting by Stellar Mass 

The most massive galaxies will produce detectable 
lensing effects over a larger area of the sky. Thus we 
choose one of the weighting methods to be M n where 
M is the object's stellar mass and n is some positive in- 
teger. The rationale for this scaling is that at the high 
mass end of the galaxy stellar mass function, the relation 
between stellar mass and total mass is non-linear. Ac- 
cording to, e.g., weak lensing, clustering, satellites, and 
abundance matching studies, the total mass increases 
with stellar mass faster than linear in for the most mas- 



sive galaxies (Mandelbaum et al.ll2006t 


Wake et al.H2011t 


Behroozi et al.ll2010t iLeauthaud et al.l 


20121: iMore et al.l 


20111). This is consistent with the fact that the central 



sider both Mf and ^£iLT M i ■ The former will 

be denoted as Wm™ while the latter as the root sum of 
the squares Wm k ■ In this paper we explore the cases 
of n = 1,2, and 3. In this section we assume to know 
precisely the correct masses, as given by the Millennium 
catalog; in Scction[5]wc allow each mass to depend on its 
respective photometric redshift to assess the reliability of 
a M n weighting method. 

4.1.4. Weighting by Luminosity 

In practice, inferring the stellar mass of an object with 
limited observational data can be difficult, leading to 
large uncertainties. For this purpose we will also explore 
weighting by luminosity, Wl- 

4.2. K GX t PDF with New Statistics 

We now proceed with the method outlined in Scction[3l 
requiring that ( q = 0.75, 1.25, 1.5, or 2 for all aforemen- 
tioned q. Contour plots for £ 9 versus £ ga i given in Figure 
[3] show how the various relative overdensiti es ar e related. 
Because each weighting scheme differs in W q from ev- 
ery other, we would like to ensure that keeping a consis- 
tent interval width E does not affect our relative spread. 
Weighting schemes with low W q would offer a higher per- 
centage of total lines of sight than those with large W q . 

To normalize our spread we multiply each W q by ^=-- 

Thus, we generalize Equation (fl]) to the following for new 
statistics: 



W q - C^gal 



< E. 



(3) 



galaxies of massive clusters and poor groups do not typ- 
ically differ in stellar mass by orders of magnitude even 



Furthermore, W q is no longer restricted to discrete in- 
teger values, but rather a continuum. Still, we expect 
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^P s ^ at under- and overdense ( q so it is necessary 
to normalize by the inverse number of LOS. We allow 
2E bins, each of length 1, from Co^Vgai — E to ( q N ga i + E 
(as previously done in Section |3| with discretely- valued 
.ZVg a i), and define A^los as the number of LOS within a 
particular bin. We then weight each K oxt value by N * 
of its respective bin when constructing the PDF. 

The variables discussed in Section |4~T1 are the principal 
contributors to gravitational lensing; however, using a 
single variable may be too basic of an approximation. We 
can expand our definition of Equation @ to allow the 
weighted sum (W q ) to be the product of characteristics 

W q w 

t=l 2=1 

where n var is the number of variables. Figurc[4](the seven 
right-most set of bars) and Table Q] show that this leads 
to greater improvements in o K when combining 1 jr with 
the established q, though in principle this can be done 
with any combination. 

In addition to generalizing W q , we can use Equation ((3|) 
to impose multiple ( q conditions. For example, we can 
require that LOS satisfy Equation ([3]) for both galaxy 
count and \/r weighting. This is a more stringent de- 
mand, as LOS must now pass two separate cuts. Note 
that our respective £ 9 (in this example Cgai) C±) values 
need not be equivalent. However, for the sake of illus- 
tration in this paper we will assume they are. We refer 
to the number of applied conditions as n con . This im- 
poses consecutive cuts that improve the quality of the 
K cxt PDF, but reduce the number of LOS. As long as 
n con is not chosen to be large enough to introduce sta- 
tistical uncertainty, we expect a sharpening of the peak 
as remaining LOS will be more relevant to the lens of 
interest. We must, however, now expand our definition 
of TVlos to incorporate combinations of every possible 
value between ±E for each of n con . Thus each LOS does 
not simply correspond to one of 2E values, but instead 
(2£)™ c °». 

We remind the reader that C? refers to the relative over- 
density of any of the aforementioned weights (e.g. galaxy 
counts, 1/r, z, L, M). Using the prescribed method we 
measure <j K for each weighting scheme. We do this by 
constructing the PDF from LOS that satisfy Equation 
^ in three variations: (1) C q ; (2) ( q and Cgai, and; (3) 
Cqi Cgah an d (i- It is worth noting that for case (2) this 
amounts to applying the same condition twice for Cgai; 
and in case (3) we have a similar redundancy for both 
Cgai and (1. 

We expect that imposing more conditions, as in cases 
(2) and (3), would lead to a smaller width of the K oxt 
PDF. Figure 2] shows er K for these latter two cases. There 
are several features here worth noting. First, changes in 
a K for C? = 0.75, 1.25, and 1.5 are relatively small when 
compared with those for C<j = 2. Thus, it is most easy to 
detect any increase or decline in o~ K for Cj = 2. Second, if 
we assume that any change in a K for each q at Co = 2 is 
indicative of the change at lower Cj (albeit on a smaller 
scale), then we can restrict the analysis to ( q = 2 to 
determine which variables constrain n ext the most. 

We also see that in all cases a K decreases when more 



conditions are imposed, as expected. Table [T] gives the 
values for the bars in Figure [4] along with Cg for case 
(1) as mentioned above. We note that for lenses with 
C 9 < 1.5 in multiple conditions, the uncertainty in re ex t is 
reduced to < 0.03, a level that is comparable to or smaller 
than the strong lens mass modeling uncertai nty in terms 
of its impa ct on the time-delay distance (e.g. ISuvu et"al"1 
120101 12013D . As current and future su rveys are expected 
to dis cover at least hundreds of lenses (jQguri fc Marshall 
I20101 ). we expect an efficient sample for cosmographic 
studies to contain lenses with relative overdensities < 
1.5. This will ensure that uncertainties due to the LOS 
structures are subdominant. 

The spread of K ex t for ( q = 2 in Table [1] is easily ex- 
plained by Figure [3] If we look at ( q — 2 for ( q ^ C ga i we 
see that a majority of its Cgai values lie at Cgai < 2. This 
causes a shift in the overall K oxt distribution, lowering 
the median and shrinking o~ K . Therefore, such a decrease 
is not the result of improving our method to find ftext? 
but the inclusion of a large number of small K oxt values 
with Cg a i < 2. This reiterates the effectiveness of impos- 
ing Cgai and other Cj conditions, as in Figure |4] and the 
remaining columns of Table [TJ 

While all weighting methods lead to a decrease in er K , 
the lowest values are for Wm" . For the K GXt PDFs 

that include the Cgai constraint, the weighting scheme 
that leads to the tightest PDF is Cgai , C m 2 , Ci = 2 with 

the corresponding PDF width as a K = 6.0423 ± 0.0006 
(see columns 3 to 6 in Table [T|) . This is a substantial 
drop from our initial finding of a K = 0.0562 ± 0.0003 for 
Crai = 2, and a K ~ 0.065 that was obtained by lSuvu et al.1 
(j2010f ) for B1608 without applying the z < 1.4 cut. 

5. TESTING THE FIDELITY OF K cxt ESTIMATES BASED 
ON PHOTOMETRIC REDSHIFTS 

Until this point our simulated galaxy catalogs from the 
Millennium Simulation have allowed us to neglect uncer- 
tainties in rcdshift and stellar mass that would normally 
arise from observations. However, getting spectroscopic 
rcdshifts for a large sample of objects - many at z > 1 
- is difficult and expensive. It is thus prudent to focus 
spectroscopic resources on the brighter objects and those 
closer to the main deflector, while using photometric red- 
shifts for the remaining objects along the line of sight. 
Because a galaxy's estimated stellar mass and luminosity 
are dependent on its redshift, their uncertainties are sen- 
sitive to errors in z. It is necessary then to estimate the 
uncertainties in K oxt associated with obtaining redshifts 
photometrically. 

5.1. Estimating photometric redshifts 

The Millennium Simulation gives magnitudes for SDSS 
u, g, r, i, z and 2MASS J, H, K bands. We use the 
Baye s ian Photometric Redshift (hereafter BPZ; iBenftezI 
l2000t iCoe et all I2006D program to calculate the photo- 
metric redshift z p hot, which is defined as the peak of the 
rcdshift PDF, for all objects. To evaluate how K ext is 
affected by the quality of Zphot, we examine three differ- 
ent band combinations. First of all, we use ugrizJHK to 
compute what we may assume to be the best approxima- 
tion to z spoc . Secondly, we use g, r, i, and K bands, in an 
effort to strike a compromise between survey speed and 
wavelength coverage, as a measure of the effectiveness of 



Improving time-delay cosmography 



7 




Fig. 4. — Values for a K for various weighting schemes with (a) successive £ ga i and ( q conditions, and (b) successive Cgal> Cq an d Cl 

conditions as measured from numerical simulations. As discussed in the text, note that cases where multiple conditions become equivalent 
(such as the leftmost bar in the top panel, corresponding to imposing two consecutive £ ga j conditions) amount to a redundancy and no 
new information is gained. The color scheme is the same as in previous plots: purple, green, red, and blue refer to £ q = 0.75, 1.25, 1.5, and 



2, respectively. The horizontal dashed lines correspond to the lowest value of a K 
find the uncertainty on re ex t is reduced by imposing more conditions. 



obtained through the weighting schemes considered. We 



our method when only a few bands are available. Lastly, 
we compute z p hot using just g, r, and i bands to evaluate 
how our method holds under the least number of bands 
from which a redshift might be computed. Typically, 
optical bands such as g, r, and i are the most readily 
available or do not require long integration times, which 
make them ideal for large-scale surveys. 



5.2. Calculating based on galaxy number counts 

with photometric data 

The errors associated with photometric redshifts are 
expected to decrease as wavelength coverage is increased. 
However, in practice obtaining ugrizJHK is observation- 
ally expensive. Thus this section serves primarily as a 
premise for the optimal strategies associated with pho- 
tometrically determining z. 

In Figure GO we plot /c ext PDFs for iVg a i with the origi- 
nal requirement of Equation (pQ) for £ g ai = 0.75, 1.25, 1.5, 
and 2 based on the spectroscopic redshifts (solid) and 
various photometric-rcdshift estimations (dashed, dot- 
dashed, dotted). With photometric- z computed using 



all 8-bands, we recover the overall shape of the /t ext 
PDF. The K oxt PDF based on griK is shown only for 
Cgai = 2, which is nearly indistinguishable from the one 
based on all bands. It is evident that fewer photomet- 
ric bands causes a shift in K oxt toward higher values. 
This is consistent with the fact that with only gri BPZ 
tends to produce photo-z with large uncertainties and 
slightly low bias. Thus, high redshift objects are incor- 
rectly assigned z < 1.4 and vice versa, but the net ex- 
change favors an increase in objects with z p hot < 1-4. 
If we assume these underestimated high-z objects are 
unlikely to be correlated with already overdense regions 
for z spcc < 1.4 (a reasonable assumption), then a uni- 
form increment, SN ga x, is accounted for along each LOS. 
Thus our LOS of previous relative overdensities go from 
Cgai^Vgai — > Cgai^Vgai+<5-^Vgai, while our new mean becomes 
= JVglT + SN eM- Multiplying Cgai by A^, we find 
that our new count based on z p hot an d satisfying (JTJ is 

Cgai(N g ai + <5iV gal ), having SN gal (C ga i - 1) more galaxies 
than the previous selection criteria. Consequently, rel- 
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Values of a K for Q q 



TABLE 1 

2 FOR VARIOUS WEIGHTING METHODS AND CONDITIONS. 



CgabCiJ— 2 CgaliCa— 2 CgallCgiCl— 2 CgaliCiiCl 



L 

z 

M 
M 2 
M 2 

" J rss 

M 3 
A/ 3 



0.0562±0.0003 
0.0534±0.0002 
0.0464±0.0002 
0.0507±0.0003 
0.0382±0.0002 
0.0319±0.0001 
0.0449±0.0002 
0.0319±0.0001 
0.0417±0.0002 



0.0181±0. 
0.0520±0. 
0.0525±0. 
0.0511±0. 
0.0492±0. 
0.0540±0. 
0.0496±0. 
0.0551±0. 



0003 
0004 
0004 
0006 
0006 
0005 
0008 
0007 



0.0478±0. 
0.0464±0. 
0.0472±0. 
0.0466±0. 
0.0508±0. 
0.0474±0. 
0.0541±0. 



0003 
0004 
0004 
0006 
0004 
0005 
0005 



0.0440±0. 
0.0455±0. 
0.0447±0. 
0.0438±0. 
0.0451±0. 
0.0436±0. 
0.0458±0. 



0005 
0004 
0006 
0009 
0007 
0011 
0008 



0.0437±0. 
0.0449±0. 
0.0430±0. 
0.0423±0. 
0.0454±0. 
0.0425±0. 
0.0464±0. 



0004 
0003 
0005 
0006 
0005 
0008 
0005 



Notes. The a K measured from numerical simulations decrease with some of the unique condition Q q = 2 (second column) because they 



incorporate a large range of f ga j values not necessarily correspondent to £ g 
Cgal = 2 (columns 3-6). 



2. This problem is fixed by always imposing initial condition 



Values of the shifts Ak^ and Aa K with £ 9 



TABLE 2 

2 FOR 2 phot COMPUTED USING A VARIETY OF BAND COMBINATIONS. 



q 


ugrizJHK 


griK 






gri 








A ^mcd 
ext. 


Acr K 


A ^mcd 
ext 


A(T B 






d 








gal 


0.0060±0.0003 


0.0023±0.0002 


0.0066±0.0005 


0.0026±0.0005 





0155±0 


0008 





0063±0 


0005 


1 

r 


0.0075±0.0003 


0.0034±0.0003 


0.0079±0.0005 


0.0025±0.0005 





0151±0 


0008 





0051±0 


0005 


L 


0.0071±0.0008 


0.0041±0.0009 


0.0075±0.0011 


0.0038±0.0019 





0145±0 


0013 





0060±0 


0019 


z 


0.0063±0.0011 


0.0039±0.0007 


0.0127±0.0017 


0.0086±0.0013 





0219±0 


0021 





0116±0 


0013 


M 


0.0054±0.0012 


0.0022±0.0011 


0.0025±0.0009 


0.0020±0.0015 





0168±0 


0011 





0065±0 


0015 


M 2 


0.0067±0.0022 


0.0029±0.0028 


0.0050±0.0008 


0.0028±0.0022 





0142±0 


0008 





0042±0 


0022 


M 2 


0.0103±0.0008 


0.0047±0.0021 


0.0105±0.0027 


0.0052±0.0023 





0202±0 


0025 





0076±0 


0023 


M 3 


0.0047±0.0015 


0.0041±0.0029 


0.0159±0.0059 


0.0085±0.0021 





0148±0 


0025 





0068±0 


0021 


M 3 

rss 


0.0080±0.0012 


0.0064±0.0015 


0.0124±0.0046 


0.0033±0.0019 





0190±0 


0012 





0085±0 


0019 


L 
r 


0.0077±0.0012 


0.0041±0.0007 


0.0077±0.0013 


0.0029±0.0016 





0115±0 


0015 





0032±0 


0016 


z 

r -> 


0.0025±0.0007 


0.0017±0.0009 


0.0047±0.0008 


0.0026±0.0007 





0143±0 


0011 





0045±0 


0007 


M 2 
T l. 


0.0061±0.0009 


0.0028±0.0007 


0.0077±0.0011 


0.0021±0.0015 





0122±0 


0015 





0043±0 


0015 


M rs B 
r i 


0.0116±0.0010 


0.0049±0.0007 


0.0121±0.0017 


0.0049±0.0017 





0203±0 


0017 





0072±0 


0017 


Mf 
<:-, 


0.0085±0.0007 


0.0047±0.0009 


0.0264±0.0037 


0.0131±0.0017 





0128±0 


0018 





0036±0 


0017 


Mf s , 
& 

r 


0.0109±0.0013 


0.0053±0.0017 


0.0100±0.0019 


0.0071±0.0021 





0136±0 


0023 





0051±0 


0021 


0.0056±0.0006 


0.0015±0.0009 


0.0017±0.0011 


0.0006±0.0017 





0194±0 


0012 





0051±0 


0017 



ative galaxy (over/under)densities based on photomet- 
ric redshifts are associated with more (over/under)dense 
lines of sight compared to spectroscopy-based galaxy 
overdensities with the same nominal value. If one were 
to ignore the difference between photometric and spec- 
troscopic redshifts, this would induce a bias in the esti- 
mation of K e xt- This bias is small when using ugrizJHK 
or even griK, but notable when only gri is available. It 
is possible, however, to further reduce such a bias by im- 
proving the photo-z to remove the small bias or by using 
exactly the same method of redshift determination in the 
actual observations and the simulations. 

The accuracy of the ^ ex t PDF should reflect the ef- 
fectiveness of a band combination's ability to estimate 
correctly an object's redshift. To quantify this accuracy 
we define the change in PDF width Acr K = crP hot - cr* ruo 
where crP hot and crj™ e are the uncertainty in K ext for z p hot 
and Zgpoc, respectively. Similarly, we define Ak™/ as the 
difference between the photometrically-determined ° d 
and the spectroscopically-determined Kcxt d - We expect 
these quantities to be the smallest for z p hot with all 8 
bands, and to increase as fewer bands are used in con- 



structing the redshift. The first row of Tablc[2]shows that 
with all 8 bands or with even only griK, the change in the 
«cxt PDF is <0.007, corresponding to <0.7% impact on 
the time-delay distance. We thus conclude that the min- 
imal set of filters necessary to achieve a 1% precision and 
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Fig. 5. — Kext PDFs constructed from LOS satisfying Equation 
GJ for Cgal = 0.75 (purple), 1.25 (green), 1.5 (red), and 2 (blue) for 
spectroscopic (solid) and photometric redshifts as measured from 
numerical simulations. Photometric redshifts are computed using 
all available optical and infrared bands (dashed), griK bands (dot- 
dashed; shown for f~ a i = 2 only) and gri bands (dotted). 
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accuracy on time-delay distance (see \Suvu et al$i20T$i, for 
a summary of cosmological implications) is three optical 
filters plus K. 

5.3. Impact of photometric redshifts on accuracy and 
precision of K cxt estimates using multiple weights 

Because Section [4] demonstrated that using character- 
istic features of galaxies provides a sharper PDF we need 
to explore how photometric redshifts affect these weights 
and their respective PDFs. Specifically, we would like to 
confirm our intuition that these weighted PDFs behave 
in the same manner as their spectroscopic counterparts, 
so that we may choose a universal optimal weighting 
method that is independent of how an object's attributes 
are obtained. 

We use our photometric redshifts for W z and rescale 
the Millennium Simulation masses and luminosities by 

( dL(z Phot j ) wriCrc dh (z) is the luminosity distance at 
rcdshift z. This is a reasonable approximation for stellar 
mass, as to first order it scales proportional to luminos- 
ity. Table [2] lists the changes in the median AK™ d = 

^ed.phot _ C ed, S pec ^ sprcad ^ = <wt <pcc q£ 

the convergence distributions. We find that for nearly 
all q, Aer K increases with fewer bands, which is consis- 
tent with our observations for £ ga i. However, when using 
either all bands or griK, the changes are below 0.01 for 
a majority of the q. 

We thus conclude that, as in the previous section, 
one should use as much information as possible to infer 
•P(ftext) f° r the observed line of sight. As shown in Fig- 
ure [3l observables like position, luminosity, redshift, and 
stellar mass add valuable information and can improve 
both the precision (by reducing a K ) and the accuracy 
(by shifting K™° d closer to the "true" value) of the infer- 
ence. In case spectroscopic redshifts are not available, we 
recommend using at least three optical bands and one in- 
frared band for the weighting schemes considered in this 
paper. For the current level of cosmological precision and 
accuracy, using griK is sufficient to constrain re ext almost 
as well as with spectroscopic redshifts. 

6. ILLUSTRATING THE METHOD WITH THE CASE 
STUDY B1608+656 

The previous sections outlined a new approach for de- 
termining K cxt for an arbitrary lens given that sufficient 
characteristics are known to calculate the relevant C, q . 
In this section we illustrate how the method works in 
practice using B1608+656 as a case study. The data 
on the B1608+656 field includes deep HST imaging in 
F555W and F814W (9 and 11 orbits, respectively; GO- 
10158, PI Fassnacht), as well as more shallow imaging 
in Gunn g, r, and i obtained with the Palomar 60- 
Inch Tel escope (full details of the observations can be 
found in lFassnacht et aT1l2006f ). This section should be 
taken as an illustration only, since the data in hand for 
B1608+656 are not sufficient to achieve the full poten- 
tial of our method. Therefore we do not revisit the 
cosmological implications of B1608+656 in this work. 
Work is in progress to collect the necessary photome- 
try and spectroscopy and future papers will present im- 
proved estimates of K GX t and cosmological parameters for 
B1608+656 and other systems. 



0.040 




Fig. 6. — Kext PDFs for £ ga i and £i values as measured from 
B 1608+656. 



6.1. Field Preparation and W q 

As a reference, we use the central portion (1.1 x 1.1 
deg 2 ) of the COSM OS field (Cosmic Evolution Survey; 
iScoville et al.| [2007t) as a sample for measuring the aver- 
age number of galaxies and also the average properties 
of the features for all lines of sight, i.e., W q . COSMOS 
data, like B1608+656 data, has ACS F814W photom- 
etry that is sufficiently deep to satisfy the upper limit 
of the I magnitude cut . We use the 2006 ACS Catalog 
(|Leauthaud et al. 2007) and matc h the galaxies to th ose 
in the 2006 Photometry Catalog (jCapak et al.l 120071) in 
order to obtain redshifts and stellar masses. We con- 
sider two objects in opposing catalogs to be identical if 
they have angular separation < 0.5". Of the ~ 124, 000 
objects in the ACS databank with 18.5 < tofsi4 < 
24.5, approximately 118,000 or 95% have photometric 
coordinate-matched counterparts. Only ~ 60 ACS ob- 
jects are found to be < 0.5" to two different objects from 
the photometry catalog. In these few cases of ambiguous 
identification, th e stellar mass from the multiband pho- 
tometric catalog (|Ilbert et al.ll2010T) is assigned to objects 
in the ACS catalog proportional to their fractional con- 
tribution to the total flux. 

Because a small subsample of the ACS catalog does not 
have matches in the photometric database, we expect our 
mean relative overdensity values for all but £ ga i and £i/ r 
to be slightly underestimated. As a solution, a correction 
factor b = N AG s/N phot is applied where A^acs and N phot 
are the total and matched number of galaxies in the ACS 
catalog, respectively. Thus multiplying the average num- 
ber of galaxies that are found in the Photometry Catalog 
by b results in the true average. This is not so apparent 
with other weighting methods, as simply multiplying W q 
by the inverse fraction of galaxies detected with ground- 
based photometry assumes the missing subset is repre- 
sentative of the entire ACS catalog. Nonetheless, given 
that the Nacs and iVphot only differ at the 5% level, we 
expect the effects of such an assumption to be small. 

To determine W q for B1608+656, we identify ob- 
j ects in the field of B1608+656 using SExtractor 
(jBertin fc Arnoutsl [l996h . To ensure a fair compar- 
ison with COSMOS, we use only a single HST or- 
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TABLE 3 

Statistics of B1608+656 relative to COSMOS 

~q ~ q W q £ kS£? ~* 

(COSMOS) (B1608+656) 

gal 41.6 83 1.997 0.0808 0.0562 2123167 

i 3.1799 6.2928 1.979 0.0741 0.0471 527734 

z 15.6181 30.6207 1.961 0.0784 0.0534 328798 

| 1.1947 2.3658 1.980 0.0695 0.0464 314931 



Notes. Column 1 lists the statistic, Column 2 are the average weights from lines of sight in COSMOS, Column 3 are the weights for 
B1608+656, and Column 4 are the relative overdensities (the ratio of Column 3 to Column 2). Columns 5 and 6 are, respectively, re e xt and 
tr K values found with all Millennium Simulation LOS that satisfy Equation for both £ ga i and £j conditions imposed. The number of 
LOS found from the Millennium data is also provided in Column 7 to indicate the reliability of the data. 



bit of imaging (so that the COSMOS and B1608+656 
images have simi lar depths), and follow the reduction 
steps outlined in iLeauthaud et ail (|2007T ). with the ex- 
ception of masking out asteroid trails, oversaturated 
stars, etc. since these are not present within 45" of 
B1608+656. Once a catalog of all objects with r < 45" 
and 18.5 < mp8i4 < 24.5 has been created, the ob- 
jects are coordinate-matched with a deeper catalog, con- 
structed from the full 11 HST orbits of imaging. This 
second catalog serves two purposes: (1) to exclude fake 
detections from the single orbit SExtractor data set, and 
(2) to produce more reliable redshifts and stellar masses 
necessary to find Q q for the detected objects in the single- 
orbit catalog. A limited number of objects (~ 15) have 
spectroscopic redshifts, or Gunn g, r, and i bands from 
ground-based imaging. Most, however, have just F814W 
and F606W photometries that make it difficult to obtain 
accurate stellar masses and photometric redshifts. 

6.2. Finding ( q for B1608 

We impose z < 1.4 for COSMOS and B1608+656. 
We find that the observed average number of galaxies 
in COSMOS is A^ = 41.6 while B1608+656 has 83 
galaxies within 45", giv ing Cgai = 1-99- This i s close to 
Cgai = 2.18 as found by iFassnacht et all (|2011l ) by com- 
paring the galaxy counts to those in pure-parallel fields. 
The slight difference in ( ga \ could be due to (1) our im- 
po sed redshift restr i ction z < z aoutce that was not applied 
bv lFassnacht"eTaLl dMl . and (2) the COSMOS field be- 
ing slightly overdense. However, it is not clear how much 
COSMOS is overdense when the redshift condition is im- 
posed. For simplicity, we neglect this correction, though 
in the future this will need to be measured before apply- 
ing it to time-delay systems for cosmological inferences. 
In a method analogous to Section 2] we compute ( q for 
each characteristic, the results of which can be seen in Ta- 
ble [3] for characteristics that are computed with a higher 
degree of accuracy (e.g., 1/r). Unfortunately the present 
data are not sufficient to estimate reliable stellar masses, 
luminosities, or accurate photometric redshifts. Deeper 
optical and NIR imaging of the field are necessary to ob- 
tain more accurate redshifts (see Tabled, luminosities, 
and stellar masses for computing £ q . 

Next we select from the Millennium Simulation lines 
of sight with the new £g values to find ft ex t- In keep- 
ing with Section 2] we impose Equation (J3J) for both £ ga i 
and £ q and find K ex t and o~ K for the resulting distribu- 
tion. These, along with the number of lines of sight, arc 
given in the last three columns of Table [3] Distribu- 
tions with Cg a i ~ ( q are closely correlated and therefore 



have large Alos- In Figure we show the K ex t PDF for 
B 1608+656 with (1) only £ ga i imposed, and (2) both £ g ai 
and Ci/r conditions imposed. The additional £i/ r condi- 
tion sharpens the PDF, leading to a decrease in a K from 
0.056 by ~ 0.01. In this specific case of B1608+656, the 
new K cxt PDF does not decrease the uncertainty on the 
final time-delay distance measurement appreciably since 
the stellar kinematics of the lens galaxy provides substan- 
tial constraints on K oxt already similar to the level that 
is achieved with the multiple £ g conditions. Nonetheless, 
the PDF of the time-delay distance is shifted to lower val- 
ues by 1 — 2% due to the lower K™ t d . We thus conclude 
that even though these effects are smaller than current 
uncertainties for a single lens, they will become impor- 
tant for reaching the ultimate goal of 1% precision and 
accuracy. 

To generalize, without the velocity dispersion as a con- 
straint, the new K oxt would have decreased the uncer- 
tainty on the resulting Hq from B1608+656 for various 
uniform cosmological priors by ~ lkms~ 1 Mpc _1 from 
~ 4kms~ 1 Mpc~ 1 . Therefore, for lens systems in which 
the lens velocity dispersion is difficult to obtain (due to, 
e.g., bright lensed images that are near the lens galaxies), 
or for very large samples for which stellar velocity dis- 
persions might not be practical, our techniques for tight- 
ening K oxt are especially valuable since the reduction in 
the uncertainty on K cxt would then translate directly to 
that on the time-delay distance (e.g., 0.01 in K ext is ap- 
proximately 1% on the time-delay distance). 

7. SUMMARY AND CONCLUSIONS 

With the goal of finding ways to measure the effects of 
the distribution of mass along the line of sight to gravi- 
tational lensing time delays, we have performed a com- 
prehensive analysis of simulated lines of sight catalogs. 
These lines of sight catalogs are based on the Millennium 
Simulation, used to compute the external convergence 
K ext via ray-tracing, as well as on semianalytic models of 
galaxy formation, used to assign observable properties to 
halos along the line of sight. 

Our main results can be summarized as follows 

1. The observed relative abundance of galaxies within 
a given aperture ( ga \ provides an estimate of re ext 
that is accurate to a few percent, depending on 
the actual under/overdensity of the observed line 
of sight. This is consistent with previous work 
(ISuvu et alJl20loh . 

2. Adding information from other obscrvables like 
stellar mass, luminosity, redshift, position of the 
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galaxies in the vicinity ol the main deflector, re- 
duces significantly the uncertainty in re GXt . The 
most significant drop in uncertainty is obtained by 
weighting each galaxy with the inverse of the pro- 
jected distance to the main deflector, followed by 
powers of the stellar mass. With this kind of infor- 
mation the uncertainty on time-delay distance aris- 
ing from K oxt can be reduced to ~ 4% from ~ 6% 
using only galaxy counts for a very overdense line 
of sight and to ~ 3% for typical lines of sight. 

3. Even though spectroscopic redshifts are valuable, 
especially for the galaxies most closely associated 
with the main deflector, photometry in three opti- 
cal bands (e.g., gri) and the near infrared (K) are 
sufficient for obtaining photometric redshifts such 
that the median and width of the K cxt change by 
< 0.007, i.e., < 0.7% on the time-delay distance. 

4. As a practical illustration, we apply this method 
to the field of B1608+656 and show that some gain 
can be made even with existing data. Better multi- 
band photometry is needed to fully realize the gains 
promised by our method. 

From these results, we conclude that with sufficient 
imaging and spectroscopy data the effects of the mass 
distribution along the line of sight on gravitational time- 
delay distances can be accounted for and the associated 
uncertainties reduced for all lines of sight. These im- 
provements - in combination with recent advances in 
the d erivation of gravitational time delays (jTewes et al.l 
I2012t ) and in the modeling of the mass distribution of 
the main deflecto r and objects in close proximity to it 
(|Suvu et all 120111) - bring us closer to the goal of 1% 
precision in cosmological distances, necessary to address 
fundamental issu es such as the nature of dark energy 
(|Suvu et al.H2012A . 

In the next decade, upcoming surveys are expected 
to deliver thousands of gravitationally lensed quasars 
(jQguri fc Marshall|[20Tol) . a number more than sufficient 
to meet the 1% goal provided effort is made to keep sys- 
tematic uncertainties under control. This will have to 
include theoretical uncertainties related to the choice of 
numerical simulations and associated semianalytic mod- 
els. Our choice of using overdensities with respect to 
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random fields, as opposed to absolute densities, mini- 
mizes the impact of the choice of this specific model. 
However, as the number of lenses with measured time 
delays increases thus reducing the observational errors 
toward the 1% level, it will be important to repeat and 
extend this study with independent cosmological simula- 
tions and galaxy formation models. This is left for future 
work. 

From an observational point of view, the future abun- 
dance of targets will change dramatically the situation 
with respect to the present time when the precision of 
time-delay cosmology is limited by the number of known 
strongly lensed quasars, and allow us to choose the tar- 
gets that give more cosmological information at fixed ob- 
servational resources. This work suggests that focusing 
follow-up efforts on specific lines of sight - those that 
are not too overdense with respect to the average of the 
universe - should result in substantial gains in efficiency. 

We would like to thank Dan Coe for his help with 
BPZ, and the anonymous referee for his/her construc- 
tive comments that greatly improved the presentation 
of this work. Z.S.G., S.H.S. and T.T. gratefully ac- 
knowledge support from the Packard Foundation in the 
form of a Packard Research Fellowship to T.T. and from 
the National Science Foundation grant AST-0642621. 
S.H. and R.D.B. acknowledge support by the National 
Science Foundation (NSF) grant number AST-0807458. 
T.E.C. is supported by an STFC studentship. P.J.M. ac- 
knowledges support from the Royal Society in the form 
of a research fellowship. C.D.F. acknowledges support 
from NSF-AST-0909119. L.V.E.K. is supported in part 
by a NWO-VIDI program subsidy (639.042.505). Sup- 
port from HST program GO- 12889 is gratefully acknowl- 
edged. This work is based in part on observations 
made with the NASA/ESA Hubble Space Telescope, ob- 
tained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS 5- 
26555. These observations are associated with program 
#GO-10158 and #GO-9822. The Millennium Simula- 
tion databases used in this paper and the web applica- 
tion providing online access to them were constructed as 
part of the activities of the German Astrophysical Vir- 
tual Observatory. 



Capak, P., Aussel, H., Ajiki, M., McCracken, H. J., Mobasher, B., 
Scoville, N., Shopbell, P., Taniguchi, Y., Thompson, D., 
Tribiano, S., Sasaki, S., Blain, A. W., Brusa, M., Carilli, C, 
Comastri, A., Carollo, C. M., Cassata, P., Colbert, J., Ellis, 
R. S., Elvis, M., Giavalisco, M., Green, W., Guzzo, L., 
Hasinger, G., Ilbert, O., Impey, C, Jahnke, K., Kartaltepe, J., 
Kneib, J. -P., Koda, J., Koekemoer, A., Komiyama, Y., 
Leauthaud, A., Le Fevre, O., Lilly, S., Liu, C, Massey, R., 
Miyazaki, S., Murayama, T., Nagao, T., Peacock, J. A., Pickles, 
A., Porciani, C, Renzini, A., Rhodes, J., Rich, M., Salvato, M., 
Sanders, D. B., Scarlata, C, Schiminovich, D., Schinnerer, E., 
Scodeggio, M., Sheth, K., Shioya, Y., Tasca, L. A. M., Taylor, 
J. E., Yan, L., & Zamorani, G. 2007, ApJS, 172, 99 

Coe, D., Bem'tez, N., Sanchez, S. F., Jce, M., Bouwcns, R., & 
Ford, H. 2006, AJ, 132, 926 

Coe, D. & Moustakas, L. A. 2009, ApJ, 706, 45 



12 



Greene et al. 



Courbin, F., Chantry, V., Rcvaz, Y., Sluse, D., Faurc, C, Tewes, 
M., Eulaers, E., Koleva, M., Asfandiyarov, I., Dye, S., Magain, 
P., van Winckel, H., Coles, J., Sana, P., Ibrahimov, M., & 
Meylan, G. 2011, A&A, 536, A53 

Courbin, P., Saha, P., & Schechter, P. L. 2002, in Lecture Notes 
in Physics, Berlin Springer Verlag, Vol. 608, Gravitational 
Lcnsing: An Astrophysical Tool, ed. F. Courbin & D. Minniti, 1 

Das, S. & Linder, E. V. 2012, ArXiv e-prints (1207.1105) 

Dobke, B. M., King, L. J., Fassnacht, C. D., & Auger, M. W. 
2009, MNRAS, 397, 311 

Ellis, R., Santos, M. R., Kncib, J.-R, & Kuijken, K. 2001, ApJ, 
560, L119 

Falco, E. E., Gorenstein, M. V., & Shapiro, I. I. 1985, ApJ, 289, 
LI 

Fassnacht, C. D., Gal, R. R., Lubin, L. M., McKean, J. P., 

Squires, G. K., & Readhead, A. C. S. 2006, ApJ, 642, 30 
Fassnacht, C. D., Koopmans, L. V. E., & Wong, K. C. 2011, 

MNRAS, 410, 2167 
Fassnacht, C. D., Womble, D. S., Ncugcbauer, G., Browne, 

I. W. A., Readhead, A. C. S., Matthews, K., & Pearson, T. J. 

1996, ApJ, 460, L103+ 
Fassnacht, C. D., Xanthopoulos, E., Koopmans, L. V. E., & 

Rusin, D. 2002, ApJ, 581, 823 
Fohlmeister, J., Kochanek, C. S., Falco, E. E., Morgan, C. W., & 

Wambsganss, J. 2008, ApJ, 676, 761 
Gavazzi, R., Treu, T., Rhodes, J. D., Koopmans, L. V. E., 

Bolton, A. S., Buries, S., Massey, R. J., & Moustakas, L. A. 

2007, ApJ, 667, 176 
Guo, Q., White, S., Boylan-Kolchin, M., De Lucia, G., 

Kauffmann, G., Lemson, G., Li, C, Springel, V., & Weinmann, 

S. 2011, MNRAS, 413, 101 
Hilbert, S., Hartlap, J., White, S. D. M., & Schneider, P. 2009, 

A&A, 499, 31 

Hilbert, S., White, S. D. M., Hartlap, J., & Schneider, P. 2007, 

MNRAS, 382, 121 
— . 2008, MNRAS, 386, 1845 

Hoekstra, H., Hartlap, J., Hilbert, S., & van Uitert, E. 2011, 

MNRAS, 412, 2095 
Ilbert, O., Salvato, M., Le Floc'h, E., Aussel, H., Capak, P., 
McCracken, H. J., Mobasher, B., Kartaltepe, J., Scoville, N., 
Sanders, D. B., Arnouts, S., Bundy, K., Cassata, P., Kneib, 
J. -P., Koekemoer, A., Le Fevre, O., Lilly, S., Surace, J., 
Taniguchi, Y., Tasca, L., Thompson, D., Tresse, L., Zamojski, 
M., Zamorani, G., & Zucca, E. 2010, ApJ, 709, 644 
Keeton, C. R. 2010, General Relativity and Gravitation, 42, 2151 
Keeton, C. R., Kochanek, C. S., & Seljak, U. 1997, ApJ, 482, 604 
Keeton, C. R. & Zabludoff, A. I. 2004, ApJ, 612, 660 
Kneib, J. P., Mellier, Y., Fort, B., & Mathez, G. 1993, A&A, 273, 
367 

Kneib, J.-P. et al. 2003, ApJ, 598, 804 

Kneib, J.-P. & Natarajan, P. 2011, A&A Rev., 19, 47 

Kochanek, C. S. 1991, ApJ, 373, 354 

— . 1996, ApJ, 466, 638 

Koopmans, L. V. E., Treu, T., Bolton, A. S., Buries, S., & 

Moustakas, L. A. 2006, ApJ, 649, 599 
Koopmans, L. V. E., & Treu, T. 2003, ApJ, 583, 606 
Lagattuta, D. J., Fassnacht, C. D., Auger, M. W., Marshall, P. J., 

Bradac, M., Treu, T., Gavazzi, R., Schrabback, T., Faure, C, 

& Anguita, T. 2010, ApJ, 716, 1579 
Leauthaud, A., Massey, R., Kneib, J.-P., Rhodes, J., Johnston, 

D. E., Capak, P., Heymans, C, Ellis, R. S., Koekemoer, A. M., 

Le Fevre, O., Mellier, Y., Refregier, A., Robin, A. C, Scoville, 

N., Tasca, L., Taylor, J. E., & Van Waerbeke, L. 2007, ApJS, 

172, 219 

Leauthaud, A., Tinker, J., Bundy, K., Behroozi, P. S., Massey, 
R., Rhodes, J., George, M. R., Kneib, J.-P., Benson, A., 
Wechsler, R. H., Busha, M. T., Capak, P., Cortes, M., Ilbert, 
O., Koekemoer, A. M., Le Fevre, O., Lilly, S., McCracken, 

H. J., Salvato, M., Schrabback, T., Scoville, N., Smith, T., & 
Taylor, J. E. 2012, ApJ, 744, 159 

Lemson, G. & Virgo Consortium, t. 2006, ArXiv Astrophysics 
e-prints 

Limousin, M., Richard, J., Julio, E., Kneib, J.-P., Fort, B., 
Soucail, G., Eh'asdottir, A., Natarajan, P., Ellis, R. S., Smail, 

I. , Czoske, O., Smith, G. P., Hudclot, P., Bardcau, S., Ebeling, 
H., Egami, E., & Knudsen, K. K. 2007, ApJ, 668, 643 



Linder, E. V. 2011, Phys. Rev. D, 84, 123529 

Mandclbaum, R., Seljak, U., Kauffmann, G., Hirata, C. M., & 

Brinkmann, J. 2006, MNRAS, 368, 715 
Meylan, G., Jetzer, P., North, P., Schneider, P., Kochanek, C. S., 

& Wambsganss, J., eds. 2006, Gravitational Lcnsing: Strong, 

Weak and Micro 
Momcheva, I., Williams, K., Keeton, C, & Zabludoff, A. 2006, 

ApJ, 641, 169 

More, S., van den Bosch, F. C, Cacciato, M., Skibba, R., Mo, 

H. J., & Yang, X. 2011, MNRAS, 410, 210 
Oguri, M. 2007, ApJ, 660, 1 

Oguri, M. & Marshall, P. J. 2010, MNRAS, 405, 2579 
Paraficz, D., Hjorth, J., & Eh'asdottir, A. 2009, A&A, 499, 395 
Pettini, M., Steidel, C. C, Adelberger, K. L., Dickinson, M., & 

Giavalisco, M. 2000, ApJ, 528, 96 
Richard, J., Stark, D. P., Ellis, R. S., George, M. R., Egami, E., 

Kneib, J.-P., & Smith, G. P. 2008, ApJ, 685, 705 
Refsdal, S. 1964, MNRAS, 128, 307 
Saha, P. 2000, AJ, 120, 1654 

Sand, D. J., Treu, T., Ellis, R. S., Smith, G. P., & Kneib, J.-P. 

2008, ApJ, 674, 711 
Schechter, P. L., Bailyn, C. D., Barr, R., et al. 1997, ApJ, 475, 

L85 

Schneider, P. 1997, MNRAS, 292, 673 
Schneider, P. & Seitz, C. 1995, A&A, 294, 411 
Scoville, N., Aussel, H., Brusa, M., Capak, P., Carollo, C. M., 
Elvis, M., Giavalisco, M., Guzzo, L., Hasinger, G., Impey, C, 
Kneib, J.-P., LeFevre, O., Lilly, S. J., Mobasher, B., Rcnzini, 
A., Rich, R. M., Sanders, D. B., Schinnerer, E.. Schminovich, 
D., Shopbell, P., Taniguchi, Y., & Tyson, N. d'. 2007, ApJS, 
172, 1 

Seljak, U. 1994, ApJ, 436, 509 

Springel, V., White, S. D. M., Jenkins, A., Frenk, C. S., Yoshida, 
N., Gao, L., Navarro, J., Thacker, R., Croton, D., Helly, J., 
Peacock, J. A., Cole, S., Thomas, P., Couchman, H., Evrard, 

A. , Colberg, J., & Pearce, F. 2005, Nature, 435, 629 
Sonnenfeld, A., Bertin, G., & Lombardi, M. 2011, A&A, 532, A37 
Sonnenfeld, A., Treu, T., Gavazzi, R., et al. 2012, ApJ, 752, 163 
Suyu, S. H., Auger, M. W., Hilbert, S., Marshall, P. J., Tewes, 

M., Treu, T., Fassnacht, C. D., Koopmans, L. V. E., Sluse, D., 
Blandford, R. D., Courbin, F., & Meylan, G. 2013, ApJ, 766, 70 

Suyu, S. H., Marshall, P. J., Auger, M. W., Hilbert, S., 

Blandford, R. D., Koopmans, L. V. E., Fassnacht, C. D., & 
Treu, T. 2010, ApJ, 711, 201 

Suyu, S. H., Treu, T., Blandford, R. D., Freedman, W. L., 
Hilbert, S., Blake, C, Braatz, J., Courbin, F., Dunkley, J., 
Grecnhill, L., Humphreys, E., Jha, S., Kirshner, R., Lo, K. Y., 
Macri, L., Madore, B. F., Marshall, P. J., Meylan, G., Mould, 
J., Reid, B., Reid, M., Riess, A., Schlegel, D., Scowcroft, V., & 
Verde, L. 2012, ArXiv e-prints (1202.4459) 

Tewes, M., Courbin, F., Meylan, C, Kochanek, C. S., Eulaers, E., 
Cantale, N., Mosquera, A. M., Magain, P., Van Winckel, H., 
Sluse, D., Cataldi, G., Voros, D., & Dye, S. 2012, ArXiv 
e-prints (1208.6009) 

Treu, T. 2010, ARA&A, 48, 87 

Treu, T., Gavazzi, R., Gorecki, A., Marshall, P. J., Koopmans, 
L. V. E., Bolton, A. S., Moustakas, L. A., & Buries, S. 2009, 
ApJ, 690, 670 

Treu, T., & Koopmans, L. V. E. 2002, MNRAS, 337, L6 
Treu, T. & Koopmans, L. V. E. 2002, ApJ, 575, 87 
Wake, D. A., Whitaker, K. E., Labbe, I., van Dokkum, P. G., 
Franx, M., Quadri, R., Brammer, G., Kriek, M., Lundgren, 

B. F., Marchesini, D., & Muzzin, A. 2011, ApJ, 728, 46 
Williams, K. A., Momcheva, I., Keeton, C. R., Zabludoff, A. I., & 

Lehar, J. 2008, ApJ, 672, 733 
Wong, K. C, Keeton, C. R., Williams, K. A., Momcheva, I. G., & 

Zabludoff, A. I. 2011, ApJ, 726, 84 
Wucknitz, O. 2002, MNRAS, 332, 951 



